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Abstract. 1RXS J 170849-4009 10 is one of four Anomalous X-ray Pulsars which emit persistent luminous radiation in soft 
X-rays (<10 keV) as well as in hard X-rays (>10 keV). In this paper we present detailed spectral and temporal characteristics 
over the whole X-ray band. For this purpose data have been used from INTEGRAL, RXTE and XMM-Newton. The hard X-ray 
(> 10 keV) time-averaged total spectrum, accumulated over four years with the imager IBIS-ISGRI onboard INTEGRAL adding 
up to 5.2 Ms net exposure, can be described by a power law with a photon index T = 1.13 ± 0.06 and extends to ~ 175 keV. The 
20-175 keV flux is (7.76 ± 0.34) x 10"" erg cirrV which exceeds the 2-10 keV (unabsorbed) flux by a factor of -2.3. No 
evidence for a spectral break is found below 300 keV. Also, no significant long-term time variability has been detected above 
20 keV on time scales of 1 and 0.5 year. Pulsed emission is measured with INTEGRAL up to 270 keV, i.e. to much higher 
energies than the total emission, with a detection significance of 12.3cr (20-270 keV). The pulse profiles from 0.5 keV up to 
270 keV show drastic morphology changes below ~20 keV. Three different pulse components can be recognized in these pulse 
profiles: 1) a hard pulse peaking around phase 0.8 which contributes to the pulse profiles above ~4 keV, 2) a softer pulse which 
peaks around phase 0.4 not contributing in the hard X-ray domain and 3) a very soft pulse component below 2 keV. A combined 
time-averaged pulsed spectrum (2.8-270 keV) from INTEGRAL, RXTE-PCA and HEXTE (collected over nine years) can be 
described with a soft and a hard power-law component: T s = 2.79 ± 0.07 and Th = 0.86 ± 0.16. In the pulsed spectrum extracted 
from a 25.5 ks net exposure XMM-Newton observation we find a discontinuity between 2 keV and 3 keV. Above these energies 
the spectrum is consistent with the spectrum taken with RXTE-PCA. The pulse profiles and the total-pulsed spectrum prove to 
be stable over the whole nine-years time span over which the data have been taken. Also detailed phase-resolved spectroscopy 
of the pulsed emission confirms the long-term stability as the spectra taken at different epochs connect smoothly. The phase- 
resolved spectra reveal complex spectral shapes which do not follow the shape of the total-pulsed spectrum. The spectral shape 
gradually changes with phase from a soft single power law to a complex multi-component shape and then to a hard single power 
law. The spectrum switches from a very hard (r = 0.99 ± 0.05) to a very soft (r = 3.58 ± 0.34) single power-law shape within 
a 0.1 -wide phase interval. The discontinuity measured between 2 keV and 3 keV with XMM-Newton is a result of a curved 
component. This component which is most apparent within phase interval 0.7-0.9 significantly contributes in the energy range 

between 4 keV and 20 keV. It has a very steep spectrum below 5 keV with a photon index T 1.5. From the phase-resolved 

spectra we identify three independent components with different spectral shapes which together can accurately describe all 
phase-resolved spectra (2.8-270 keV). The three shapes are a soft power law (T = 3.54), a hard power law (T = 0.99) and 
a curved shape (described with two logparabolic functions). The phase distributions of the normalizations of these spectral 
components form three decoupled pulse profiles. The soft component peaks around phase 0.4 while the other two components 
peak around phase 0.8. The width of the curved component (~0.25 in phase) is about half the width of the hard component. After 
4U 0142+61, 1RXS J1708-40 is the second anomalous X-ray pulsar for which such detailed phase-resolved spectroscopy has 
been performed. These results give important constraints showing that three dimensional modeling covering both the geometry 
and different production processes is required to explain our findings. 
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In the past five years, several new developments have changed 
our view on Anom alous X -ray Pulsars (AXPs) drastically (see 
for a recent review Kaspi 2007). AXPs are young isolated ro- 
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tating neutron stars. Their spin periods and spin derivatives 
indicate that the surface magnetic fields of AXPs range from 
~6 x 10 13 G to ~7 x 10 14 G (assuming a dipole field; B = 
3.2 x 10 19 VpPG). The X-ray luminosity of persistent AXPs 
can be orders of magnitude higher than explained by spin-down 
power. Nowadays, it is generally accepted that AXPs are pow- 
ered by ma gnetic dissipation and that AXPs belong to the mag- 
netar class (Duncan & Thompson! 1992 |: Thompson & Duncan 



1991 1 19951 119961: iThompson et al.l l2002l) . To the magne 
tar class also belong Soft Gamma-ray Repeaters (SGRs) 
dWoods & Thompson! BOOli iMereghetti etal] 12007b . So far, 
there are in total 13 confirmed magnetars known. 

Traditionally, AXPs are soft X-ray sources (~0.5-10keV) 
with thermal-like spectra (kT ~ 0.4 - 0.7 keV) plus a 
soft power-law-like component with a photon index F ~ 
3-4. This traditional view changed since the detections 
of extremely bright and persistent h ard X-ray emission 
from AXPs by INTEGRAL Q20 keV iMolkov et al.l 12004 ; 
Revnivtsev et all 12004 Iden Hartog et all |2004 . The hard X- 



ray spectra of these AXPs can be described with power-law 
models with pho ton indices T ~ 0.9 -1.4 and are observed 
up to 230 keV dKuiper et alj l2004at Iden Hartog et alj [2006; 



Kuip er et alj|2006t Iden Hartog et alj |2008). Because no coun 
terparts have been found in the MeV gamma-ray domain by 
CGRO-COMPTEL (0.75-30 MeV), t hese spectra are expected 
to show breaks. For AXP 4U 0142+6l lde"n Hartog et all 00081) 
measured the first spectral break in an AXP hard X-ray spec- 
trum. The maximum luminosity of 4U 0142+61 is estimated to 
be near 280 keV. 

Moreover, iKuiper et all d2004al 120061) found that a sig- 
nificant part of the hard X-ray emission is pulsed with 
spectral shapes harder than the total-emission spectra. For 
4U 0142+61 pulsed emissi on was detected up to 160 keV 
by Iden Hartog et al.l (120081) . These authors performed for 
4U 0142+61 a phase-resolved spectral analysis of the pulsed 
emission (i.e. without the dominating non-pulsating DC emis- 
sion) in a broad band from 0.5 keV to 160keV exploiting data 
from four space missions. In that studv lden Hartog et alj Q2008) 
showed that in the pulse profiles at least three different pulse 
components can be recognized which appeared to have differ- 
ent spectra. Therefore, these pulse components are genuinely 
different from each other. These findings provide new con- 
straints for modelling the underlying physical processes re- 
sponsible for these spectra and the geometry of the produc- 
tion sites in the AXP magne t ospheres. Theoretical models cur- 
rently in development (e . g. iBeloborodov & Thompson! 2007 : 
iHeyl & Hernquistll2005bl: iBaring & Harding! l2007b ~do so far 
not explain these detailed results. 

In this paper we present new high-energy results from 
AXP 1RXS J 170849-4009 10 (hereafter 1RXS J1708-40). 
1RXS J 1708-40 has a spin period of P = 1 1 .0 s and a p eriod 
derivative of P = 1.9 x 10~ n ss" 1 dGavriil & Kaspill2002l) . The 
characteristic age (t = P/2P) dipolar surface magnetic-field 
strength of 1RXS J1708-40 inferred from these values are 9.0 
kyr and 4.7 x 10 14 G, respectively. In the 2-10 keV band its 
luminosity Lx is 1 x 10 35 erg s " 1 , similar to all other AXPs as 



Soon after its discovery dSugizaki et alJI 19971) 1RXS J1708-40 
has been included in th e RXTE monitoring program of AXPs 
( Gavriil & Kaspil 2002 ). providing a long base-line of obser- 
vations with a wealth of new information li ke timing parame- 
ters, pulsed-flux history and glitches (see e.g. lKaspi et aljl999t 
Dibetal] l2008). For this source an intensity-hardness correla 



tion in the soft X-ray band (<10 keV) is claimed (IRea et al 



2007 : Campana et al.l2007h . An indication th at this correlation 



also holds in the hard X-ray domain (>10 keV lGotz et al . 2007) 
is not confirmed in this work. 

1RXS J 17 08-40 was first de t ected i n the hard X-ray regime 
(>20 keV) bv lRevnivtsev et alJ d2004 in a deep INTEGRAL 



Galactic centre map. [Kuiper et al. (2006) discovered pulsed 
hard X-ray emission from this source using RXTE-PCA 
and HEXTE data. These authors also analysed early-mission 
INTEGRAL data adding up to almost 1 Ms effective on-source 
exposure and found that the hard X-ray spectrum could be de- 
scribed by a power-law model with a photon index of 1 .44 + 
0.45. Timing analysis revealed significant pulsed emission in 
the 20-300 keV band. 

This work is a co ntinuation of the approach followed in 
den Hartog etail d2008l) for AXP 4U 0142+61 applying the 



suming a distance of 3.8 kpc (Du rant & van Kerkwiikll20063) . 



same methods now for 1RXS J1708-40. We address again the 
total high-energy window, using data collected over four years 
of INTEGRAL observations adding up to ~5 Ms of on-source 
exposure, nine years of RXTE monitoring adding up to an ex- 
posure of ~600 ks and a ~25 ks XMM-Newton observation. A 
significantly improved time-averaged INTEGRAL total spec- 
trum is presented, as well as those for shorter (~ 1 year and ~0.5 
year) data sets to look for possible long-term time variability. 
Most importantly, we present the (time-averaged) total-pulsed 
spectra as well as (time-averaged) phase-resolved pulsed spec- 
tra from INTEGRAL, RXTE and XMM-Newton. 



2. Observations and analysis 

2.1. INTEGRAL 

The INTErnational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) is a hard X-ray/soft gamma-ray (3 keV - 
8 MeV) mission which has been operational since October 
2002. The payload consists of two main coded-mask imaging 
telescopes. 1) I BIS, Imager on Board the INTEGRAL Satellite 
dUbertini et al.ll2003b . has a wide field of view (FOV) of 29° x 
29° (full- width zero response) and it has an ~12' angular res- 
olution. The low-energy dete ctor is called the IN TEGRAL Soft 
Gamma-Ray Imager (ISGRI: ILebrun et"aL 2003 ) and it is sen- 
sitive between ~20 k eV and ~300 keV 2) The Spectrometer 
for INTEGRAL (SPI; IVedrenne et ai1l2003l) has a better spec- 
tral resolution but only a moderate imaging capability with an 
angular resolution of ~2°5. SPI is sensitive between ~20 keV 
and ~8 MeV and has a FOV of 35° x 35° (full-width zero re- 
sponse). 

For the best sky map reconst ruction IN TEGRAL performs 
observations in dither patterns dJensen et alJl2003l) besides the 
occasional staring observations. Typical pointings (Science 
Windows, ScWs) last -1800-3600 s. Due to the limited vis- 
ibility windows on the Galactic-centre region, the observations 
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Table 1. Summary of the INTEGRAL observations of 1RXS J1708-40. The following information is given: the revolution 
intervals, the observations time spans both in MJD as in calendar dates, the number of ScWs, exposure times (t exp ), and the 
effective on-source exposure (f e //). 



Rev. 








Time 


span 






ScWs 


^exp (ks) 


ten (ks) 


037- 


485 


52668 


-54013 


Jan. 


29, 


2003 


- Oct. 5, 


2006 


5161 


11959.1 


5191.0 


037- 


063 


52668 


-52751 


Jan. 


29, 


2003 


- Apr. 22, 


2003 


706 


1394.3 


700.7 


100- 


120 


52859 


- 52920 


Aug 


8, 


2003 


- Oct. 8, 


2003 


717 


1843.6 


709.0 


164- 


185 


53051 


-53115 


Feb. 


16, 


2004 


- Apr. 20, 


2004 


592 


1317.0 


519.0 


222- 


246 


53225 


- 53298 


Aug 


8, 


2004 


- Oct. 20, 


2004 


698 


1590.3 


618.7 


283- 


307 


53407 


-53481 


Feb. 


6, 


2005 


-Apr. 21, 


2005 


1027 


2256.9 


1116.4 


345- 


366 


53594 


- 53655 


Aug 


12, 


2005 


- Oct. 12, 


2005 


470 


1009.6 


450.5 


399- 


431 


53756 


-53850 


Jan. 


21, 


2006 


- Apr. 25, 


2006 


517 


1405.8 


534.0 


466- 


485 


53960 


-54013 


Aug 


13, 


2006 


- Oct. 5, 


2006 


434 


1141.6 


542.6 



form sets of consecutive orbital revolutions (Revs). In Table Q] 
a list of these sets are given with the total exposure time for the 
selected ScWs for which 1RXS J1708-40 was within an an- 
gle of 14°.5 from the pointing direction of IBIS, as well as the 
effective on-source exposure, which is reduced due to off-axis 
viewing angles. 

The region where 1RXS J1708-40 is located is often ob- 
served with INTEGRAL as it is ~15° from the Galactic centre. 
After four years of observations the total exposure on this re- 
gion adds up to 12 Ms taken in 5161 pointings, after screening 
for Solar flares and erratic count-rates due to the passages of 
the spacecraft through the Earth's radiation belts. Taking into 
account the off-axis angles the effective on-source exposure is 
5.2 Ms. 

2.1 .1 . INTEGRAL spectra from spatial analysis 

The spectral analysis of the IBIS-ISG RI da t a follo ws the same 
procedures as applied in Iden Hartog et al.l (120080 . The shad- 
owgrams of 5161 ScWs are deconvolved into sky images in 
20 energy intervals using the Off-line Scientific Analysis soft - 



ware package version 5.1 (OSA; see Gol dwurm et al.l 12003). 
The energy intervals are exponentially binned between 20 keV 
and 300 keV. At the source position of 1RXS J 1708-40 aver- 
aged count rates are determined for all energy bins. Photon 
fluxes are generated by normalizing the AXP count-rates to 
the Crab count-rates determined from INTEGRAL Crab obser- 
vations during Rev. 102, and the Crab spectrum. The spectral 
shape of the Crab spectrum we adopted is a curved power-law 
shape 



1.5703(14)x(£ r /0.06335) 



-2.097(2)-0.0082(16)xln(£ y /0.06335) 



(1) 



Table 3. Summary of the RXTE observations of 1RXS J1708- 
40. The following information is given: the observation IDs, 
time spans of the observations in calendar dates, PCA exposure 
times for PCU-2 in units ks. 



Obs. ID. 




Time 


span 




?PCA 


ALL 


Jan. 12, 


1998- 


Nov. 16, 


2006 


606.512 


Set A 


Jan. 12, 


1998- 


Oct. 26, 


2003 


310.480 


SetB 


Nov. 1, 


2003 - 


Nov. 16, 


2006 


296.032 






Set A 








30125 


Jan. 12, 


1998- 


Jan. 08, 


1999 


59.896 


40083 


Feb. 6, 


1999- 


Mar. 11, 


2000 


52.568 


50082 


Apr. 21, 


2000- 


May 12, 


2001 


34.576 


60069 


May 6, 


2001 - 


Feb. 20, 


2002 


24.544 


60412 


May 20, 


2001 - 


May. 23, 


2001 


9.928 


70094 


Apr. 2, 


2002- 


Mar. 20, 


2003 


55.600 


80098 


Apr. 16, 


2003- 


Oct. 26, 


2003 


73.368 






SetB 








80098 


Nov. 1, 


2003- 


Feb. 23, 


2004 


26.416 


90076 


Feb. 28, 


2004- 


Feb. 28, 


2005 


126.784 


91070 


Mar. 6, 


2005 - 


Feb. 25, 


2006 


79.000 


92006 


Mar. 5, 


2006- 


Nov. 16, 


2006 


63.832 



been extracted for every energy bin using a maximum likeli- 
hood fi tting procedure th at considers the instrumental response 
of SPI (lKn6dlsederll2004l) . Also here the measured flux is nor- 
malized to the Crab spectrum, which has been obtained using 
SPI observations of Revs. 43-45. 

2.1 .2. INTEGRAL timing analyses 

For the IBIS-ISGRI timing analysis of A XP 1RXS J 1708-40 



(Eq. 2 in iKuiper et aD l2006h . where F y is expressed in we followed ear l ier ap plied procedures dKuiper et al. 2006; 

•,2 



ph/(cm s MeV) and E y in MeV. For fitting these (and other) 
spectr a the software package X-spec version 12.3 (lArnaud 
19961) and private software are used. All errors quoted in this 
paper are lcr errors. 

For the INTEGRAL observations up to Rev. 366 SPI spec- 
tral analysis is performed above 140 keV, from where the sen- 
sitivity of SPI is comparable or higher than ISGRI. In the 
140 keV - 1 MeV energy range the data were binned in 4 en- 
ergy bins. From the total-exposure sky maps source fluxes have 



den Hartog et al .1 |2008). All available data taken between 
Rev. 037 and 485 were used (see Table [TJ. Photons were se- 
lected only from non-noisy detector elements that could be il- 
luminated by the source through the open mask elements by 
more than 25%. After correcting for ins trumental, onboard pro- 
cessing and ground-station time delays ( Walter et alj2 003), the 
event arrival times are Solar-system barycentered (JPL DE200 
Solar- system ephemeris) adopting the best known source po- 
sition dlsrael et alj|2003l) . The barycentered events are folded 
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Table 2. Phase-coherent ephemerides - derived from RXTE-PCA monitoring data - valid for the analysed INTEGRAL observa- 
tions. 



AXP 


Start 


End 


INTEGRAL 


*. Epoch 


V 


V 


V 


O 




[MJD] 


[MJD] 


range (Revs) 


[MJD, TDB] 


[Hz] 


xlO- 13 [Hzs- 1 ] 


xl0~ 22 [Hz s- 2 ] 




1RXS J 1708-40 


52590 


52960 


009-133 


52590 


0.0908982198657 


-1.58268 


0.00 


0.4333 


1RXS J 1708-40 


53050 


53325 


163-255 


53050 


0.0908919724436 


-1.62994 


2.75 


0.5025 


1RXS J1708-40 


53321 


53520 


254-320 


53377 


0.0908875116829 


-1.49563 


-34.1 


0.7559 


1RXS J 1708-40 


53555 


53690 


332-377 


53555 


0.0908852012859 


-1.77561 


36.4 


0.8587 


1RXS J 1708-40 


53741 


54055 


395-499 


53741 


0.0908825979045 


-1.57228 


-0.813 


0.1842 



using a phase-connected ephemeris (see Table [2] and Sect. 12.21 
for the ephemeris creation). The (TDB) time to pulse phase 
conversion taking into account consistent phase alignment for 
each ephemeris is provided by the following formula: 

1 ,1 , 

O(f) = v-(f-f Epoch ) + -v-(f-f Epoch ) / + -v-(f-f Epoch ) J -<Do.(2) 

For the pulse profiles, detection significances are estimated 
by applying a Z^ tes t on bin-free pulse-phase distributions. 
(iBuccheri et al.lfl983l) . Truncated Fourier series with two har- 
monics are used to fit the pulse profiles, where the minimum of 
these fits are defined as the DC level. 

The total-pulsed spectrum is derived by obtaining the num- 
ber of excess counts above the DC level from each pulse profile 
for all considered energy bands. These excess counts are nor- 
malized in Crab units (by applying the same procedure for the 
Crab pulsar) and are converted i nto photon fluxes us ing a Crab- 
pulsar spectrum model (Eq. 3 in lKuiper et al.l :2006): 



F y = 0.4693(21)x(£ r /0.04844) 



■1.955(7)-0.071 0(78)xln(£ r /0.04844) 



(3) 



The analysis to create phase-resolved spectra (Sect. 13.3.21 ) 
is similar. However, for the phase-resolved spectra only the ex- 
cess counts within selected phase intervals are used. 

2.2. RXTE 

For the timing analysis of the INTEGRAL data the regular 
monitoring observations with the Rossi X-ray Timing Explorer 
(RXTE) du ring the INTEGRAL observations are of great 
importance (iGavriil & Kaspil 120021) . because the INTEGRAL 
source count rate is too low to independently detect the pul- 
sation. This is possible with the data fr om the Proportional 
Counter Array (PCA; Ijahoda et al.1ll996h aboard RXTE. The 
PCA is a non-imaging instrument sensitive in the 2-60 keV en- 
ergy range and it consists of five collimated units (PCUs 0^4) 
with a ~1° FOV. Using these RXTE-PCA monitoring data we 
created accurate phase-connected ephemerides valid during the 
INTEGRAL observations. The details of the ephemerides can 
be found in Table|2] The pulsar phase was set arbitrarily to zero 
where the pulse profile for energies between ~2.5 and 10 keV 
reached a sharp mini mum (see e.g . Fig .|7]B). This is a different 
choice than made in iKuiper et al] (120061) . resulting in a phase 
shift in the profiles of ~0.2. 

For the study of the pulsed emission of 1RXS J 1708-40 
with RXTE-PCA, we have used all publicly available data 
spanning nine years of observations. The data are analysed in 



two sets. Se t A co ntains the data which have been used by 
Kuiper et al. (2006) and were taken over a period of almost 
six years. Most of these data were taken before the start of the 
INTEGRAL observations. Set B covers the remaining part of 
the publicly available data taken during INTEGRAL operations 
and span about three years. All data used in this work have been 
listed in Table [3] It can be seen that Set A and B have similar 
total exposures on 1RXS J 1708-40. 

The procedu res for the RXTE spectral-timing analysis 
as outlined by Kuiper et alJ d2006l) . and also applied in 
den Hartog et al. (2008), have been followed. The pulse pro- 
files for any selected energy band can be described suffi- 
ciently accurate by truncated Fourier series with three harmon- 
ics above constant DC levels. The excess (pulsed) counts above 
these DC levels can then be converted into flux units using PCU 
exposure-weighted response matrices (see e.g. Sect. 3.2 of 
Kuiper et al.ll2006l) . To derive unabsorbed PCA time-averaged 
total-pulsed and phase-resolved spectra, a Galactic absorption 
column jVh = 1.47 x 10 22 cirT 2 was used (fixed value), derived 
from XMM-Newton spectral fits (see Sect. 13. 1 .3l >. 

The High-Energy X-Ray Tim ing Experiment (HEXTE; 15- 
250 keV; iRothschild etaL1ll998b aboard RXTE has also been 
used to extend the PCA spectral-timing study to higher en- 
ergies. HEXTE consists of two independent detector clus- 
ters (0 and 1) each with a -1° FOV. HEXTE allows for de- 
tailed consistency checks with the INTEGRAL results ob- 
tained for similar energy windows. Due to the co-alignment 
of HEXTE and the PCA, the same long monitoring observa- 
tions of 1 RXS J 1 708 -40 can be used for both instruments . Fully 
consistent procedures have been applied for the spectral-timing 
analysis, including the use of the ephemerides listed in Table[2] 
and following the same procedures for deriving pulsed counts 
and fluxes of the pulsed emission. On 1RXS J1708-40 the total 
deadtime corrected exposures for clusters and 1 are 239.04 
ks and 259.50 ks, respectively. 

2.3. XMM-Newton 



XMM-Newton (Ja nsen et al.H2001l) has been operational since 
early 2000. Onboard are three CCD cameras for X-ray imaging, 
namely two E PIC (European Photon Imaging Ca mera) MPS 



dTurner et al. 2001 ) cameras and one EPIC-PN dStriider et al 
200 ll) camera. All cameras have a FOV of ~30' and are sensi- 
tive in the energy range ~0.3-12 keV. 

We have (re)analysed a publicly available 45 ks XMM- 
Newton observation on 1RXS J 1708-40. The data were taken in 
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Table 4. Summary of the power-law model fits to the 
INTEGRAL-IBIS-ISGRI total-spectra of 1RXS J1708-40 
summing observations taken in intervals of ~1 year and for the 
total exposure. For each spectrum are given the revolution in- 
terval over which 1RXS J1708-40 observations are summed, 
the photon index, the 20-150 keV flux and^ . 



Rev. 


r 




F20-150 X 10- 1J 


X\ (d0f) 








[erg cm" 2 s 4 ] 




037-120 


1.18 ± 


0.14 


5.47 ± 0.43 


1.64 (14) 


164-246 


1.40 ± 


0.15 


5.97 ± 0.55 


1.09 (14) 


283-366 


0.97 ± 


0.12 


7.13 ± 0.44 


0.52(14) 


399^185 


1.00 ± 


0.21 


8.72 + 0.61 


1.67 (14) 


037^185 


1.13 ± 


0.06 


6.61 ± 0.23 + 


1.13 (14) 


^20-175 = 


(7.76 ± 


0.34) x 10-" erg cm" 2 s 


-i 



the night of August 28 to 29, 2003 (Obs. Id 0148690101). The 
data are analysed using SAS v. 7.0 and the latest calibration 
files that were available (May 2007). Next, we considered the 
EPIC-PN which operated in small-window mode with medium 
filter. In this mode the maximum time resolution is 6 ms. These 
data are checked for solar (soft proton) flares by creating a light 
curve with events with energies larger than 10 keV. A count- 
rate distribution was created from this light curve to which a 
Gaussian was fitted to determine the mean count rate of the 
high-energy photons. Good-Time Intervals (GTIs) were created 
allowing only time stamps for which the count rate was lower 
than 0.078 counts per second which corresponds to the fitted 
mean count rate plus three times the width of the distribution 
(3<t). The GTIs add up to 25.5 ks of exposure. All single and 
double events, i.e. patterns less than and equal four, are selected 
within the energy range 0.3 keV to 12 keV. The arrival times of 
the selected events are barycentered. 

There was no need to correct for pile-up effects as the count 
rate of 1RXS J1708-40 was low enough for the PN in small 
window mode with medium filter. To extract the source counts 
a circular extraction region with a radius of 35" was used. As 
this source is not so bright it was possible to select a back- 
ground extraction region from the small window with a radius 
of 45". No out-of-time event correction has been applied as the 
effect is approximately 1%. The extracted spectrum has been 
binned oversampling the energy resolution by a factor of three 
and then rebinned occasionally (for energies >5 keV) to ensure 
a minimum of 25 counts per bin. 

For timing analysis, the barycentered events are folded us- 
ing an appropriate ephemeris (ephemeris 1 in Tabled. To cre- 
ate pulsed spectra pulsed excess counts are extracted as de- 
scribed in Sect. 12.1.21 fitting the pulse profiles with truncated 
Fourier series with the first three harmonics above a constant 
background. 

3. Results 

In this section we present results from INTEGRAL, XMM- 
Newton and RXTE. First, we present the total (pulsed plus DC) 
spectrum using INTEGRAL data above 20 keV and look for 
possible long-term time variability. Then we derive the com- 
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Fig. 1. Total-flux (pulsed + DC) spectra of 1RXS J1708-40 col- 
lected during four one-year intervals of INTEGRAL observa- 
tions, corresponding to revolution intervals (from left to right); 
037-120, 164-246, 283-366 and 399-485. The upper limits are 
2cr. The model lines are single power-law fits (see Table|4j>. 

plementary XMM-Newton total spectrum below 12 keV. RXTE 
is a non-imaging mission and due to the presence of nearby 
stron g sources and large gradi ents in the Galactic ridge emis- 
sion Jvalinia & Marshalll l998). it can in this work only be used 
for studying the timing signal of 1RXS J1708-40. The pulse 
profiles are then presented and compared for all three missions. 
Finally, the results from spectral timing analyses in the form of 
total-pulsed spectra and phase-resolved spectra for all data sets 
are presented and discussed in detail. This leads to the identi- 
fication of genuinely different components contributing to the 
pulsed high-energy emission. 

3.1. Total spectra and long-term variability 

3.1 .1 . INTEGRAL ISGRI persistent time-averaged 
spectra and long-term variability 

Listed in Table Q] are eight sets of INTEGRAL observations, 
each covering about half a year. 1RXS J1708-40 is not bright in 
hard X-rays, therefore we first performed the spectral analysis 
averaging observations over one-year time intervals to obtain 
significant flux values till above 100 keV. The resulting four 
spectra are shown in Fig.Q] They can be described satisfacto- 
rily with single power-law models. The fit results are listed in 
TableE] The power-law indices range from 0.97 ±0. 12 to 1 .40± 
0. 1 5 and the 20-1 50 keV fluxes range from (5 .47 ± 0.43) x 1 0" 1 1 
to (8.72 ± 0.61) x 10~ n erg cm^s -1 which can be seen as a in- 
dication for time variability. When we compare these yearly 
averages with the long-term four year average (see Sect. 13.1.21 , 
we find that the spectra taken in Revs. 037-120 and 399-485 
are most deviating from this time-averaged spectrum, i.e. at the 
~2cr and ~3cr level, respectively (see error contours in Fig.|2|. 

Alternatively, the standard deviations (s = 
yjl/(n - 1) Yji( x i ~ x) 2 ) f° r me power-law indices and the 
20-150 keV fluxes relative to the weighted mean are 0.20 and 
1.46 x 10~ n erg cm _2 s _1 , respectively. Therefore, the power- 
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Fig. 2. Error contours (lcr, 2<x and 3<x) of the power-law model 
fits to the total spectrum (20-300 keV) for observations 037- 
120 and 399-485 compared to the total time-averaged fit (037- 
485). Out of the four observations (each averaged over one 
year) these two differ most from the four-year time-averaged 
spectrum. 



law shape is stable within 18% (lcr) and the 20-150 keV flux 
is stable within 22% (lcr). 

To explore the indication for variability further, we also fit- 
ted the spectra of the half-year long time sets (Table [TJ plus a 
more recent short, lower-statistics observation (Revs 534-542, 
Feb. 27 - Mar. 25, 2007, 173 ScWs, 307 ks effective on-source 
exposure) with single power-law models, but over the narrower 
energy band 20-70 keV. The results are plotted in Fig. [3] The 
20-70 keV flux measurement for Revs. 399—43 1 is the only one 
that deviates at a 3<x level from the time-average value. Since 
this is one single 3cr deviation for nine measurements, we con- 
sider this only as an indication for time variability. 

The bottom panel of Fig. [3]shows no indication for si gnif- 
icant vari a tion in spectral index, either. Ilsrael et aL ( 2007t) and 
Dib et al.l (120081) reported the detection of glitches and candi- 
date glitches of 1RXS J1708-40 during the INTEGRAL oper- 
ations, the epochs are indicated in Fig. [3] This figure shows 
that there is no indication for an obvious change in spec- 
tral slope nor flux triggered by one of the glitches. Possibly 
there is a change in spectral slope after the glitch on MJD 
53366. However, the difference between T = 1.27 + 0.10 and 
r = 0.98 ± 0.08, the average values before and after the glitch, 
is not significant (2.3cr). 



Gotz et alJ d2007l) analysed part of the INTEGRAL data 



used in this work, they did not use Revs. 164-185, 345-366 
and 399-431, and claim the detection of long-term hard X- 
ray (20-70 keV) variability based on an indication for lower 
flux values in the INTEGRAL measurements after the sec- 
ond candidate glitch (epoch marked in Fig. [3]), correlated with 
lower fluxes measured with Swift in the 1-10 keV band (their 
Fig. 1). However, our more complete analysis including full 
spectral fitting of the data in stead of using merely count rates 
does not confirm the reported trend in the INTEGRAL data for 
1RXS J1708-40 (see Fig.©. Therefore, we do not support the 



: i 



o 

O t— 
CM v 



+- 



H 1 1 1 1 H 



5.3x1CT 



5.35x1Cf 
MJD 



5.4x1CT 



Fig. 3. Fluxes (20-70 keV) and photon indices for the half-year 
time sets specified in Table Q] and one shorter observation (see 
Sect. I3.1.U . Indicated as horizontal dotted lines are the time- 
averaged values wi th their lcr errors ( grey bands). Gli tches and 
candidate glitches (Ilsrael et alJl2007tlDibet al.ll2008h are indi- 
cated as vertical black and grey lines, respectively. 



claimed correlated variability between the hard X-ray and soft 
X-ray emission from 1RXS J 1708-40. 



3.1 .2. INTEGRAL time-averaged total spectrum 

Combining all the data listed in Table Q] a time-averaged 
spectrum is extracted using 5.2 Ms on-source exposure time 
(Fig. g). 1RXS J 1708-40 has been detected with high signif- 
icance up to 175 keV with a 4.2cr detection significance in 
the 152-175 keV band. Fitting this spectrum with a single 
power law yields an excellent fit (xj. = 1.13 with 14 dof) with 
a photon index F = 1.13 ± 0.06 and a 20-150 keV flux of 
(6.61 + 0.23) x 10- 11 erg cirrV 1 . Its 20-175 keV luminosity 
is L?,o-n.s = 1.34 x 10 35 ergs" 1 assuming a distance of 3.8 kpc 
(Dur ant & van Kerkwiik 2006a). It is remarkable that there is 



no hint for a spectral break up to the highest energies in our 
INTEGRAL-ISGRI spectrum, particularly when compared to 
the r eported COM PTEL upper limits for energies above 750 
keV dKuiper et al.l l2006), which impose a spectral break be- 
tween 100 keV and 750 keV (see Fig. g). 

Unfortunately, it turned out that 1RXS J1708-40 is 
not bright enough above ~100 keV to be detected with 
INTEGRAL-SPI, even for this long exposure. In Fig. g] three 
2cr upper limits derived from SPI spatial analysis are added. 
These limits are all above the extrapolation of the power-law fit 
to the ISGRI flux values, contrary to the case of 4U 0142+61, 
for which SPI uppe r limits provided evide nce for the presence 
of a spectral break (Iden Hartog et al ] l2008h . 
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Fig. 4. High-energy spectra of 1RXS J1708-40. In this figure 
the following is plotted: the unabsorbed total spectra of XMM- 
Newton (<12 keV) and INTEGRAL (with triangl e markers) in 
black, also in black three COMPTEL upper limits dKuiper et al 
l2006h . three INTEGRAL-SPI upper limits in grey (with trian- 
gle markers); also in grey a power-law fit to the INTEGRAL- 
IBIS spectrum, in blue a logparabolic fit to the INTEGRAL- 
IBIS, SPI and COMPTEL data, total pulsed spectra of XMM- 
Newton in black, RXTE-PCA and HEXTE are shown in blue 
and aqua, and the total pulsed spectrum of INTEGRAL-ISGRI 
in red (with triangle markers). 



Following Iden Hartog et al.l d2008l) we have fitted all 
INTEGRAL and COMPTEL spectral information (including 
limits) with a logparabolic function; 



(4) 



where Eq (in units keV) is the pivot energy to minimize cor- 
relations between the parameters and Fq is the flux (in units 
ph cm _2 s _1 keV _1 ) at Eq. This function is a power-law if the 
curvature parameter j3 is equal to zero. Assuming this spectral 
shape we get an acceptable broad-band (20 keV - 30 MeV) fit 
with best-fit parameters a = 1.637 ± 0.049,jS = 0.261 ± 0.035 
andFo = (1.68±0.08)xlO- 6 phcm- 2 s- 1 keV- 1 at£ = 143.276 
keV. The peak energy £p ea k is 287^75 keV. This value lies re- 
markably cl ose to the peak energy found for 4U 0142+61 (i.e. 
279^1 keV: lden Hartog et alJEooi) . In Fig.Hboth the power- 
law and the logparabolic fit are drawn. 

3.1 .3. XMM-Newton total spectrum 

For energies below 12 keV we extracted the absorbed to- 
tal (pulsed + DC) spectrum using XMM-Newton EPIC-PN 
data (see Sect. 12.3b . In order to obtain an estimate for the 
Galactic absorption column (Nh) we fitted the spectrum glob- 
ally with a canonical logparabolic function, including fixed 
INTEGRAL parameters for the hard X-ray contribution above 
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Fig. 5. IBIS-ISGRI pulse profile of 1RXS J 1708-40 between 20 
keV and 270 keV. This profile has a 12.3cr significance using 
a Z 2 test (fit shown as a solid curve). The fitted DC level is 
indicated with an horizontal line. The grey lines and the colours 
indicate three phase intervals Phi, II and III (see Table |6j. The 
colours are consistently used in this paper in figures showing 
results of phase-resolved analyses (see Sect. 13 .3 .21 . 



~8 keV. We derive an N u of (1.47 + 0.02) x 10 22 cm" 2 , 
which can be com pared with the v alue (1.36 + 0.03) x 10 22 
cm obtained by [Rea et al. d2005h fitting the same XMM- 
Newton data with an absorbed black-body plus a power- 
Si 



law model. iDurant & van Kerkwiik (2006b) used a model 
independent approach analysing X-ray gratin g spectra taken 



with the Reflection Grating Spectrometer (den Herd er et al 



2001b onboard XMM-Newton. Their value for N H of (1.40 ± 
0.4) x 10 22 crrT 2 is consistent with both estimates. We adopted 
N H = 1.47 x 10 22 ctrT 2 in this work for the XMM-Newton and 
RXTE analyses. The total unabsorbed spectrum is shown in 
Fig.ffl The 2-10 keV unabsorbed flux is (3.398 ±0.012) x 10" 1 1 
erg cm _2 s _1 . The error is statistical only. The 2-10 keV unab- 
sorbed fluxes for N H = 1.40 x 10 22 and 1.36 x 10 22 cm 2 are 
(3.361 +0.009)x 10- 11 and (3.339 + 0.013)x lO" 11 erg cnrV 1 , 
respectively. These values are within 2% of our value. 

3.2. Pulse profiles 

3.2.1 . INTEGRAL and XMM-Newton pulse profiles 

Kuiper et al. I d2006h showed for the first time pulsed hard X- 



ray emission (>10 keV) from 1RXS J1708-40 using data from 
RXTE-PCA, RXTE-HEXTE and INTEGRAL-ISGRI. For the 
INTEGRAL pulse profiles ~IA Ms on-source exposure was 
used, resulting in a 5.9<x detection for energies 20-300 keV. 
In this work, we present INTEGRAL pulse profiles using 
~5.2 Ms on-source exposure. The result is a very much im- 
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Fig. 6. Pulse profiles from soft to hard X-rays. XMM-Newton 
pulse profiles in the energy range 0.5-12.0 keV are shown in 
panels A, B, C and D. Panels E, F, G and H show INTEGRAL 
pulse profiles in the energy range 20-270 keV. The differential 
energy ranges are indicated in the figures. The phase intervals 
are indicated as in Fig. [5] 
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Fig. 7. Pulse profiles from RXTE. RXTE-PCA pulse profiles 
in the energy range 2.2-63.9 keV are shown in panels A, B, C, 
D and E. RXTE-HEXTE pulse profiles are shown in panels F, 
G and H within the energy range 14.8-267.6 keV. The phase 
intervals are indicated as in Fig. [5] 



prov ed pulse profile with a 12.3<x detection significance (Z^ 



test; IBuccherietal.ll 1983b for energies 20-270 keV (Fig. 01. 



The profile shows a single pulse which peaks around phase 0.8 
with a steep trailing wing dropping off to the DC level at phase 
~1.05 (0.05). On the leading wing there appears to be a weak 
pulse or shoulder. 

Presented in Fig.|6]E-H are four exponentially-binned dif- 
ferential INTEGRAL pulse profiles. The profiles have signifi- 
cances of 6.3cr, 6.8cr, 6.2cr and 3.5cr, respectively, using a 
test. All profiles show single pulses which all peak at phase 
~0.8, like for the total pulse in Fig. [5] 



The complementary high-statistics XMM-Newton pulse 
profiles at lower energies (0.5-12 keV) are shown in Fig.|6]A- 
D. The pulse profile with the lowest statistics, in panel D, still 
has a significance of 10.7cr (Z^ test). The drastic change in mor- 
phology moving up in energy from 0.5 keV up to 270 keV is 
evident and will be investigated further. Note that in the tim- 
ing analysis the pulsed emission has been detected to higher 
energies than the total emission in the sky maps. 
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3.2.2. RXTE-PCA and RXTE-HEXTE pulse profiles 

The RXTE-PCA bridges the observational gap (12-20 keV) in 
energy between XMM-Newton and INTEGRAL, with XMM- 
Newton being sensitive down to -0.5 keV and INTEGRAL ex- 
tending the coverage to -300 keV. In addition, RXTE-HEXTE 
is sensitive over about the same energy band as INTEGRAL- 
ISGRJ. Analysis of data from these three missions allow for 
consistency checks, particularly when the observations are 
covering the same epochs and for checks on long-term vari- 
ability. RXTE observations set A (Table covers about 6 
years ( mostly before the l aunch of INTEGRAL) and was 
used by Kuip er et al.1 (12006 ). Set B overlaps in time with the 
INTEGRAL observations up to December 2006. Therefore, we 
first compared the time-averaged RXTE-PCA pulse profiles in 
5 differential energy bands for set A with those for set B. We 
did not find significant differences in these time averaged (over 
years) profile sh apes. This is in agreement with the findings by 
Dib et al ( 20081) . who compared the profile shapes for seven 
(glitch-free) intervals of the total set A plus B. We should note 
that the PCA gain changed over time, what can offer an expla- 
nation for some indicati ons for small c hanges with time in the 



profiles as discussed in lDib et al. (2008) 



In order to exploit the maximal statistics, we created pulse 
profiles using all available data collected with the PCA and 
HEXTE over nine years (see Table [3] and Sect. 12.21 1. Before 
summing the profiles obtained at different epochs we corrected 
for the above-mentioned PCA gain drift in the conversion from 
channels to energy. In Fig.Qthe RXTE-PCA (panels A-E) and 
HEXTE (panels F-H) pulse profiles are presented. The Towest- 
significance' profiles above ~15 keV (Fig.|7]D-H) have signif- 
icances of 18.5cr,4.8<x, 12. Oct, 6.4crand4.3cr, respectively. 

3.2.3. Pulse profile changes with energy 

Comparing the profile shapes derived with INTEGRAL (Fig.|6]l 
with those obtained with the PCA and HEXTE in similar en- 
ergy bands (Fig.[7]l, shows very good agreement, INTEGRAL 
achieving the best statistics above 70 keV. Above ~ 20 keV 
there does not seem to occur significant changes in morphol- 
ogy with energy. The pulse profiles are statistically similar 
which was tested usin g a combina tion of a Pearson x 2 ( see 



e.g. Sect. 11.2 of Eadie et a l.lll971l) and a Run test (see e.g. 



Sect. 11.3.1 of lEadieetalJI 19711) 



Below 20 keV, the situation is different: The XMM-Newton 
and RXTE-PCA profiles show strong morphology changes as 
a function of energy. This has been noted in a number of ear- 



lier pu blications (e.g. ISugizaki et al.lll997t llsrael et al .1 12001 



Kuip er et al.l l2006). In fact, the morphology is changing faster 
with energy than can be seen in the broad energy bins se- 
lected for Fig. [6] and Fig. [7] In addition, there are system- 
atic effects which should be realized when comparing XMM- 
Newton and PCA profiles. For example, in the measured en- 
ergy band 2.2-3.9 keV (Fig.|7]A) one can see in the PCA pro- 
file a small secondary peak at phase -0.8, which is absent in 
the corresponding XMM-Newton profile (Fig. |6]B). This can 
be explained with the significantly coarser energy resolution of 
the PCA, accepting more events with photon energy above 4 



keV in the band with measured energies below 4 keV. For the 
same reason, this secondary peak at phase 0.8 in the PCA 3.9- 
8.2 profile (Fig.|7]B) is higher than found in the corresponding 
XMM-Newton profile (Fig.|6]C), because this pulse component 
is maximal for energies above 8 keV. Therefore, a more accu- 
rate comparison can be made between the pulse characteristics 
measured with XMM-Newton and the RXTE-PCA by perform- 
ing phase-resolved spectroscopy in which the energy resolution 
is automatically taken into account (addressed in Sect. l3.3.2l i. 

Nevertheless, the XMM-Newton and RXTE-PCA profiles 
consistently show that the soft X-ray pulse has two compo- 
nents which peak at phases -0.4 and -0.8. The pulse peaking at 
phase ~0.4 is spectrally softer than the pulse peaking at phase 
-0.8. As mentioned above, the latter is not visible in Fig. [6] A 
and B. This peak starts contributing above -4 keV (Fig. |6]C 
and Fig. 0B) and is the dominating pulse component already 
above -8 keV up to and including the hard X-ray band. 

The soft pulse (peak at phase -0.4) is gradually decreasing 
in strength as a function of energy and seems to disappear en- 
tirely at hard X-rays. Finally, there is an indication for a third 
pulse component, visible in the XMM profile below 2 keV 
where the RXTE-PCA has no sensitivity. This component at 
phase -1.0 fills up the soft pulse to a saw-tooth-like profile 
(Fig.|6]A), at least, there is a sudden, significant change in pulse 
shape in comparison to the profile above 2 keV. 

3.3. Energy spectra of the pulsed emission 

For INTEGRAL-ISGRI, XMM-Newton-PN, RXTE-PCA and 
RXTE-HEXTE pulsed spectra were created by analyzing the 
pulse profiles following the procedure described in Sect. [2] In 
Fig. [5] an example is given of how the pulsed counts are ex- 
tracted as excess counts above a flat background level. In this 
case for 1RXS J1708-40, the -40% pulsed emission in the to- 
tal INTEGRAL energy band is separated from the -60% DC 
component. For INTEGRAL pulse profiles for seven exponen- 
tially binned energy intervals within 20-300 keV were created. 
For XMM-Newton, RXTE-PCA and HEXTE pulse profiles in 
20 (0.5-12 keV), 15 (2.7-32.1 keV) and 5 (14.8-267.5 keV) 
energy intervals were produced, respectively. In this section 
we first show the total-pulsed phase-averaged spectra. Then 
we show for all three missions the results of phase-resolved 
spectroscopy applied for three broad phase bins selected in the 
INTEGRAL profile in Fig. [5] Finally, we show phase-resolved 
pulsed spectra for narrow (A<f> = 0.1) phase bins for RXTE- 
PCA and INTEGRAL leading to the identification of three 
spectrally different components which contribute to the total 
pulsed emission. 



3.3.1. Total-pulsed spectra 

In Fig. |4] the time-averaged total-pulsed spectra are shown for 
all four instruments using all available exposures. The flux 
values from XMM-Newton, RXTE-PCA, RXTE-HEXTE and 
INTEGRAL are very consistent, yielding a continuous broad- 
band spectrum from 0.5 keV up to almost 300 keV. The RXTE- 
PCA total-pulsed spectrum smoothly bridges the energy gap 
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between INTEGRAL and XMM-Newton. The RXTE-HEXTE 
spectrum is in agreement with the INTEGRAL spectrum and 
also connects to the RXTE-PCA spectrum. 

For RXTE-PCA we first created spectra for sets A and B 
separately, to verify whether the time-averaged 1RXS J1708- 
40 spectra differ in the years before and during the INTEGRAL 
observations. Within the statistical errors the spectra appeared 
to be identical, in flux and shape. T his is in agreement with the 
recent findings of iDib et al.l (120081) . who also concluded that 
the pulsed emission is very stable over the whole time span 
from 1998 until 2006. Therefore, the combined spectrum with 
maximal statistics using both sets of data is shown in Fig.|4]and 
is further described below. 

Fitting INTEGRAL, RXTE-PCA and HEXTE simultane- 
ously with two power-law components yields an excellent fit 
for energies above 2.8 keV (see Table [5]). The soft power law 
has a photon index T = 2.79 + 0.07 and the hard power law has 
a photon index F = 0.86+0.16 ixl = 0.51, dof = 22). Fitting the 
INTEGRAL spectrum separately it can be described by a sin- 
gle power-law model with a photon index F = 0.98 ± 0.3 1 (x 2 r 
= 0.70, dof = 5), which is nicely within errors of the combined 
INTEGRAL-RXTE fit. 

To determine the model parameters of the 0.5-12.0 keV 
XMM-Newton spectrum the contribution of the hard power- 
law spectral component (significant above ~5 keV) were added 
with fixed parameters (as determined from the INTEGRAL- 
RXTE fit) to the model and thereby effectively subtracted. 
Fitting the 0.5-12.0 keV XMM-Newton spectrum with a 
power-law yields a rather poor fit with photon index F = 
3.08 ± 0.03 (x 2 r = 1.62, dof = 15). This photon index is slightly 
softer than found for RXTE-PCA. However, the XMM-Newton 
spectrum in Fig.|4]shows a discontinuity between 2 and 3 keV, 
just below the RXTE-PCA energy band. Phase-resolved spec- 
troscopy in the next section will show that this discontinuity 
in the total-pulsed spectrum between 2 and 3 keV is genuine. 
Indeed, for energies above 2.8 keV an excellent fit can be ob- 
tained with a power law with a photon index F = 2.87 ± 0.12, 
which is in agreement with the RXTE-PCA spectrum above 
2.8 keV (in flux and index, see Table |5). The normalization of 
the XMM-Newton spectrum below 2.8 keV is higher than the 
one above 2.8 keV, but the spectral shape is the same with an 
photon index F = 2.89 ± 0.06. 

The XMM-Newton results imply an average pulsed frac- 
tion (defined as the pulsed emission divided by the total emis- 
sion) of (35.6+1 .4)% for the 2-10 keV band. The INTEGRAL- 
RXTE fit gives a 20-150 keV pulsed flux of (2.56 + 0.40) x 
10~ n erg cm _2 s _1 , resulting in an average pulsed fraction of 
(39±6)%. The broad-band 2-10 keV and 20-150 keV pulsed 
fractions suggest a fairly constant pulsed fraction over the total 
high-energy window. However, in Fig. Hit can be seen that the 
measured pulsed flux (204-300 keV) is on the extrapolation of 
the power-law fit to the total spectrum. At these high energies 
the pulsed fraction could be as high as 100%. 



Table 5. Spectral-fit parameters of the total-pulsed 
INTEGRAL, XMM-Newton, RXTE-PCA and RXTE-HEXTE 
fits (see Sect. ED). The Nu is fixed to 1.47 x 10 22 cm" 2 . The 
integrated fluxes are in units erg cm _2 s _1 . Fq is the normaliza- 
tion at 1 keV in units ph cm _2 s _1 keV _1 . The subscripts s and h 
stand for 'soft' and 'hard' in the combined fit. 



TlNTEGRAL 

f20-150keV 

x] (dof) 

F 20-270 keV 



INTEGRAL 20-270 keV 
0.98 + 0.31 
(2.53 + 0.40) x 10-" 

0.70 (5) 
(4.89+ 1.12) x 10-" 

XMM-Newton 0.5-2.8 keV f 



r ' XMWl-Newton 
F(} XMM-Newton 

Xl (d0f) 



2.89 ± 0.06 
(1.95 + 0.07) x 10" 



0.46 (7) 

" . 5 _ 2kcV (4.69 ± 0.24) x 10-" 

XMM-Newton 2.8-12 keV 1 " 



^XMM-Newton 
FoXMM-Newton 

X 2 r (d0f) 

F 1 

1 2-10keV 



2.87 + 0.12 
(1.43 + 0.25) x 10~ 2 

0.39 (6) 
(1.21+0.05) x 10-" 



RXTE-PCA+HEXTE+INTEGRAL 2.8-270 keV 



r h 

^h,20-150keV 

X 2 r (dof) 

^2-lOkeV 
F 10-20 keV 
f20-150keV 
-f*20-270kcV 



2.79 ± 0.07 
(1.34 + 0.11) x 10~ 2 

0.86 + 0.16 
(2.40 + 0.32) x 10-" 

0.51 (22) 
(1.243 + 0.01 l)x 10-" 
(0.333 ± 0.010) x 10 " 
(2.60 + 0.35) x 10-" 
(5.16+ 1.06) x 10-" 



f : Model includes fixed parameters r h and F h 2 o-i50kev of the 
combined RXTE-INTEGRAL spectrum to correct for the 
contribution of the hard X-ray spectrum in the XMM-Newton band. 
The integrated flux is the model flux including the hard component. 

3.3.2. Phase-resolved pulsed spectra 

The morphology of the 20-270 keV INTEGRAL pulse profile 
(Fig. [5} was used to define three broad phase intervals: Phi 
contains the DC level (green coloured in figures); Phil con- 
tains the shoulder in the INTEGRAL band (blue coloured in 
figures); and Phlll contains the main pulse in the INTEGRAL 
band (magenta coloured in figures). The phase interval bound- 
aries are listed in Tableland are indicated with vertical lines 
and colour bars in Figs. |5] |6] and |7] In the latter figures is also 
visible how these phase intervals relate to the profile shapes at 
energies below 10 keV. 

In the upper panels of Fig. |8]the resulting pulsed-emission 
spectra for these three phase intervals are shown, in each case 
together with the total-pulsed spectrum, and for XMM-Newton 
and INTEGRAL also with the total-emission spectra; XMM- 
Newton in the upper-left panel; RXTE-PCA in the top-center 
panel; and in the top-right panel XMM-Newton, RXTE-PCA 
and INTEGRAL together. In the RXTE-PCA analysis we have 
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Fig. 8. Phase-resolved results for XMM-Newton and RXTE-PCA are presented in the left and center panels, respectively. These 
data and the INTEGRAL (overlapping) results are shown together in the right panels. Total (except for RXTE), total-pulsed and 
phase-resolved pulsed spectra are presented in the top panels. The total spectra and total-pulsed spectra are plotted in black. The 
phase resolved spectra for Ph I, Ph II and Ph III (see Table |6]l are plotted in green, blue and magenta, respectively. The XMM- 
Newton pulsed spectra are also indicated with open circle symbols while the INTEGRAL pulsed spectra are indicated with a 
triangle symbols. The data points with arrows indicate that the flux values have significances less than 1.5<x and the positive lcr 
error is drawn. In the bottom panels rj - defined as the fraction of the total pulsed emission within a phase interval - is shown in 
the same colour scheme as the spectra. 



Table 6. Selected phase intervals for extraction of high-energy 
spectra, using the pulse-shape morphology of the INTEGRAL- 
IBIS-ISGRI 20-270 keV pulse profile (Fig.©. 





Phase interval 


Component 


Phi 


[0.95,0.25) 


DC level INTEGRAL profile 


Phil 


[0.25,0.65) 


Shoulder in INTEGRAL band 


Phlll 


[0.65,0.95) 


Main hard INTEGRAL pulse 



again verified that the three phase-resolved spectra for Set A 
and Set B are fully consistent: there is no sign for long-term 
variability in shape and flux down to the 5%-level. 

We first focus on the pulsed-emission spectrum of Phlll 
which contains the main hard INTEGRAL pulse. This Phlll 
spectrum follows the total-pulsed spectrum in the INTEGRAL 
band, as expected. However, Ph III shows below 10 keV a com- 



pletely different spectral behaviour with respect to the total- 
pulsed spectrum. Below 10 keV where the total -pulsed spec- 
trum turns up towards lower energies, this spectrum keeps go- 
ing down steeply. The pulsed emission of Ph III reaches a mini- 
mum in luminosity around ~2.5 keV. At lower energies it turns 
up again. This remarkably spectral shape is measured consis- 
tently with XMM-Newton and RXTE-PCA, as can be seen 
in the right panel of Fig. [8] (magenta coloured data points). 
Therefore, we can exclude a systematic problem as cause of 
this appearance. It is also remarkable how smoothly the com- 
bined XMM-ftwton and nine-year-average RXTE-PCA spec- 
tra join the four-year-average INTEGRAL spectrum. This is 
strong evidence for a very stable configuration of the produc- 
tion sites of the high-energy emission in the magnetar magne- 
tosphere over almost a decade. 

It is obviously not straightforward to fit the spectrum of 
Ph III. We attempted to make an empirical description to 'quan- 
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Table 7. Spectral-fit parameters for Ph III INTEGRAL, XMM- 
Newton, RXTE-PCA and RXTE-HEXTE fits (see Sect. l332t , 
The Nu is fixed to 1.47 x 10 22 crrT 2 . kT and £break are m units 
keV. Fq is the normalization at 1 keV in units ph cnr 2 s _1 keV _1 . 
The integrated fluxed are in units erg cm _2 s _1 . 



XMM-/VewfoK+RXTE-PCA+HEXTE+ 
INTEGRAL 0.5-270 keV 



kT 

FokT 

n 

^break 

r n 
r,„ 

F20-150keV 

X\ (dof) 

F . 5 -2keV (ergcirrV 1 ) 
(ergcm-V 1 ) 
F 2 o-i50kev(ergcm" 2 s" 1 ) 
F2o-270kcv(ergcm" 2 s _1 ) 



0.23 ± 0.02 
(2.82 ± 0.64) x 10- 5 
-0.55 ±0.14 
4.92 ±0.11 
1.77 ±0.09 
(2.05 ± 0.52) x 10- 5 

0.02 ± 0.50 
(0.86 ± 0.22) x 10-" 

0.90 (29) 
(0.19 ± 0.07) x 10-" 
(0.23 ± 0.07) x 10-" 
(1.54 ± 0.18) x 10-" 
(3.74 ± 0.65) x 10-" 
INTEGRAL 20-270 keV 

0.16 



/(erg cm 2 s ') 



1.09: 

(1.53 ± 0.12) x 10-" 



xl (dof) 1.15(5) 
F 20 -270kev(erg cnrV 1 ) (2.81 ± 0.42) x 10" 11 



tify' the discontinuities using all data points from INTEGRAL, 
RXTE-PCA and XMM-Newton. We succeeded in describing 
the spectrum over the whole 0.5-270 keV band with a black 
body and three power-law segments (i.e. one broken power law 
plus a power law; = 0-90 dof = 29, see Table |7}. The total 
fit contains a soft thermal black-body component below 2 keV 
(kT = 0.23 ± 0.02 keV). This component can be due to a tail 
of the soft pulse which peaks in Ph II. Above ~2 keV the spec- 
trum drastically hardens. The spectrum up to the break energy 
£break = 4.92 + 0.11 keV is extremely hard with a photon in- 
dex L = -0.55 + 0.14. After the break a softer power law is 
required with photon index r = 1.77 ± 0.09. Finally, in order 
to account for the harder INTEGRAL spectrum a third power 
law is required with photon index F = 0.02 + 0.50, which is 
dominating above ~60 keV. It is obvious that very different 
production processes contribute in this phase interval (Ph III) 
to the high-energy (pulsed) spectrum. 

Fitting only the INTEGRAL spectrum with a power law, 
a photon index F = 1.09 + 0.16 gives the best description. 
(TableE). 

Also the spectral shape of the pulsed emission in Ph II (cen- 
tered on the soft pulse peaking at phase ~0.4, blue coloured in 
figures) is very different from that of the total-pulsed spectrum. 

The independent XMM-Newton and RXTE-PCA spectra 
are again fully consistent. A broad soft component below about 
8 keV turns in the RXTE band above 8 keV into a compo- 
nent with power-law index ~2, but is not detected around 20 
keV. The flux measurements close to energies as high as 100 
keV are likely not related to the softer X-ray components and 
can be explained as emission from a tail of the main hard-X- 




Energy [keV] 

Fig. 9. The pulsed fractions of 1RXS J1708-40 as a function 
of energy using XMM-Newton below 10 keV and INTEGRAL 
above 20 keV. The total pulsed fraction is drawn in black and 
plotted in colour are the pulsed fractions for the three phase 
intervals (see Tabled. 



ray INTEGRAL pulse in Ph III. No sensible model fits can be 
made to the total combined Ph II spectrum. A good fit can be 
obtained considering only XMM-Newton data, namely with a 
broken power-law with break energy Streak = 2.49 ± 0.15 keV, 
consistent with the energy where in the total pulsed spectrum 
the discontinuity was found. Below the break energy the power- 
law model has a photon index F = 2.83 + 0.04, above this en- 
ergy F = 3.38 ±0.06. 

Ph I contains the left wing of the soft pulse below 8 keV, 
and harbours the DC level in the INTEGRAL band (green 
coloured in figures). Therefore, no significant flux is measured 
at hard X-rays. Also in this case no sensible model combination 
(black body, logparabola, power law) can describe the com- 
bined XMM-Afewfon-RXTE-PCA spectrum. Furthermore, the 
XMM-Newton spectrum exhibits again a drastic drop in flux 
between 2 and 3 keV. The fact that the RXTE-PCA data points 
between 2 and 3 keV do not follow this sudden drop, might be 
explained with the coarser energy resolution of RXTE-PCA, 
combined with a strong gradient in the sensitive area and un- 
derlying photon spectrum. 

A different way to study the relative contributions of the 
pulsed components in the three selected phase intervals as a 
function of energy is by plotting the r\ parameter, defined as 
the fraction of the total-pulsed emission. It shows more clearly 
discontinuities in the dependence on energy. In the bottom pan- 
els of Fig. [8] 77 is presented for the different observations. The 
bottom panel for XMM-Newton shows that the Ph III contribu- 
tion to the pulsed emission is a few percent up to ~2.5 keV. 
Above 2.5 keV the contribution increases as a function of en- 
ergy to ~50% at 10 keV. This picture is fully confirmed by 
RXTE, which shows that r\ is still increasing above 10 keV up 
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Phase 

Fig. 11. Phase distributions of the normalizations Fq of the 
three spectral-model components (see Table [8) taken at pivot 
energy Eq = 15.0433 keV for the 10 narrow phase-interval 
spectra. Indicated with crosses and a thin grey line are the nor- 
malizations of the soft power-law component, with triangles 
and thick grey line the hard power-law component and with 
squares and a black line the composite spectral component. 

to ~60% around 30 keV and stabilizes around this value in the 
INTEGRAL band. 

The opposite can be seen for the contribution to the pulsed 
flux from Ph II. A high r\ of ~70%, fairly constant up to ~3 keV 
from where the contribution abruptly decreases to ~40% at 
10 keV. Also this decline is observed with RXTE. The be- 
haviour in the INTEGRAL band is not clear. Finally, Ph I shows 
a somewhat similar trend as Ph II with a significant contribu- 
tion below ~3 keV of ~28%, which decreases to ~8% above 
~8 keV and only upper limits in the INTEGRAL band. 

For XMM-Newton and INTEGRAL the pulsed fraction is 
shown in Fig. [9] The pulsed fraction (total-pulsed emission di- 
vided by total emission) of 1RXS J1708-40 is around -40% 
from soft to hard X-rays (see also Sect. l3.3.Tb . However, we can 
note here, that the pulsed fraction is not constant with energy, 
but shows significant variations with energy below 10 keV. In 
the INTEGRAL band the error bars are too large to draw firm 
conclusions. In Fig. [9] the pulsed fractions for the three phase 
intervals are also shown (i.e. emission in Phi, II or III divided 
by total emission). As was also clear from the dependence of 
j] on energy in Fig. [8] the contribution of the hard X-ray pulse 
increases with energy above ~3 keV, while the pulsed fractions 
of the other two phase intervals are decreasing with energy. 

3.3.3. Narrow-band phase resolved pulsed spectra 
from RXTE 

The phase-resolved spectroscopy presented in Sect. 13.3.21 se- 
lecting three broad phase intervals, unambiguously revealed 



that spectrally different components contribute to the total 
pulsed emission of 1RXS J1708-40. A similar conclusion wa s 
earlier reached for AXP 4U 0142+61 dden Hartog et alJl2008h . 
However, the situation appears to be much more complex for 
1RXS J 1708-40. The XMM-Newton spectra below ~2 keV 
suggest a soft thermal component, which is not unexpected. 
However, this component can not be easily separated from 
components with apparently complex spectral shapes con- 
tributing to the total pulsed emission above 2 keV. Therefore, 
we exploited the high statistics of the total ~600 ks of RXTE- 
PCA data to explore in more detail the spectral characteris- 
tics as a function of phase above 2 keV and the connection 
to the INTEGRAL band. We repeated the phase-resolved spec- 
troscopy but now for 10 narrow bins each covering 0. 1 in phase. 
Applying the method described in Sect. 12.21 the 10 phase- 
resolved spectra have been extracted and are shown in Fig.fTUl 

The variations in spectral shape with phase are drastic. In 
the first 4 phase bins between phase 0.0 and 0.4, we see a spec- 
tral tail of the soft component peaking in energy in the XMM- 
Newton band below 1 keV. These spectra can be accurately fit- 
ted for energies above 2.8 keV with a single power-law model. 
The index for the best fit to the sum of the spectra between 
phase 0.0 and 0.3 is F = 3.54 ± 0.01. In the phase bins 0.4- 
0.5 and 0.5-0.6 a significant hardening of the spectra sets in 
for energies between about 5 keV and 20 keV. In the following 
phase bins 0.6-0.7 and 0.7-0.8, getting closer to the phase of 
the main INTEGRAL pulse, evidence is visible for the pres- 
ence of the hard power-law component with index ~1.0 above 
20 keV. This latter hard power-law spectrum is dominating the 
energy budget for the last three phase bins. More surprisingly, 
in the two phase windows between 0.7 and 0.9, centered on the 
INTEGRAL main pulse, a very hard part of the spectrum is ap- 
parent in a narrow energy interval above ~2 keV. This is most 
extreme for the phase interval 0.7-0.8 with a 'local' photon in- 
dex F = -1.49 + 0.21 (notice the sign!) up to a break energy 
£Weak = 4.69 + 0.04. Finally, the spectrum for phase interval 
0.9-1 .0, the last fraction of the trailing wing of the INTEGRAL 
main pulse, can be fitted with a single power law with best-fit 
photon index F = 0.99 ±0.05. Moving only 10% in phase to the 
next phase interval 0.0-0.1, the hard spectrum 'switches' and 
becomes again very soft with T = 3.58 ± 0.34. 

From Fig.[l0]it is evident that different spectral components 
contribute to the total pulsed emission and dominate in differ- 
ent phase intervals. Two components can clearly be identified: 
A soft power-law component with index F ~ 3.6 and a hard 
power-law component with index F ~ 1.0. Then there appears 
to be a third component, dominating in the energy band ~3 to 
20 keV for phases between 0.7 and 0.9. A clear characteristic 
of this spectral component is a very steep rise above 2.3 keV 

(F 1.5). We investigated what shape such a third component 

should have in order to explain, together with two power-law 
components, all spectra in Fig. \W\ This shape should give a 
clear constraint for any proposed production mechanism. 

We followed the following procedure. The hard power-law 
component was adopted from the best fit to the spectrum for 
phase 0.9-1.0: F = 0.99. Then we investigated empirically 
what spectral shape, together with the power-law model with 
T = 0.99, can accurately describe the spectrum in phase inter- 
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Fig. 10. Narrow-band phase-resolved spectroscopy of 1RXS J1708-40 using RXTE-PCA and INTEGRAL-ISGRI data. Each 
panel shows the spectrum of a phase interval with A<fi — 0.1. The phase intervals are indicated in the panels. The data points with 
arrows indicate that the flux values have significances less than 1.5<x. The positive lcr errors are drawn. Total spectral fits and the 
three components (when contributing significantly) are drawn in black and grey, respectively (See Sect l3.3.3l and Fi g . IT0T> . 
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Table 8. Model parameters of the three spectral components 
that are used to fit the 10 narrow -band phase-resolved RXTE- 
PCA and INTEGRAL-ISGRI spectra. The three components 
are a soft and a hard power-law component, and a curved com- 
ponent which is the sum of two logparabolic functions. For 
each component Fq (in units ph cm _2 s _1 keV _I ) is the normal- 
ization taken at the pivot energy Eq = 15.0433 keV which is 
plotted as a function of phase in Fig.fTTl 



Power-law components F = Fn x 



Soft power law T, = 3.54 
Hard power law T h = 0.99 

Composite component 



FaX 



k = 7.2632 x 10~ 3 
or, = 11.635866 
/J, = 4.867669 
a 2 = 2.3099390 
p 2 = 1.6262881 



val 0.7-0.8. This turned out to be the sum of two logparabolic 
functions of which the parameters are given in Table [8] Finally, 
for the index of the soft power-law model we adopted the 
best fit to the summed spectrum in the phase interval 0.0-0.3: 
r = 3.54. 

We now attempted to describe the measured spectral distri- 
butions in the 10 narrow phase intervals in terms of just these 
three model shapes with free normalizations. Interestingly, the 
resulting total fits are very satisfactory for all narrow phase in- 
tervals (see Fig. [Toll. The normalizations of the three models are 
shown in Fig. [TTJ Effectively, this figure shows the pulse pro- 
file of each spectral component contributing to the total-pulsed 
spectrum, i.e. a complete decoupling based on the measured 
spectral and timing characteristics. In Fig.QJJwe show for all 
10 spectra the contributions of the three spectral components 
to the model fit, but only when components were statistically 
required to reproduce the measured spectra. Fig. [TTJ together 
with Fig. [TTJ show now the characteristics of the three compo- 
nents which can reproduce together the total pulsed emission, 
the pulse shapes and their spectra. Furthermore, we have shown 
above that our results indicate that these characteristics are in- 
variant over almost a decade of observations. 



4. Summary 

In this paper we have presented new and more detailed char- 
acteristics of 1RXS J1708-40 in the hard X-ray/soft gamma- 
ray regime (>20 keV) using INTEGRAL-ISGRI and SPI 
data. Followi ng the approach in ou r earlier paper on AXP 
4U 0142+61 dden Hartog et al.ll2008l) we extended the energy 
window to lower energies to obtain a broader high-energy view 
using archival RXTE and XMM-Newton observations. The em- 
phasis of this work on broad-band phase-resolved spectroscopy 
was aimed at the identification of distinctly different compo- 



nents contributing to the total-pulsed emission. What follows 
is a summary of the main results. 

4.1. Total high-energy emission above 20 keV of 
1RXS J 1708-40 

1) Using all available INTEGRAL data taken with INTEGRAL 
ISGRI during 2003-2006 - adding up to 5.2 Ms effective on- 
source exposure - the time-averaged 20-175 keV spectrum can 
be described by a power-law with photon index F = 1.13±0.06. 
The luminosity is 1.34 x 10 35 erg s~'(20-175 keV) adopting a 
distance of 3.8 kpc (Sect. 1X7721 FigM Table^. 

2) There is no indication for a spectral break in the INTEGRAL 
spectrum up to energies of ~ 175 keV. Including in the spec- 
tral fit earlier published COMPTE L flux upper limits for ener- 
gies above 750 keV dKuiper et al.l 120061) and assuming a log- 
parabolic shape for the spectrum above 20 keV, the maximum 
luminosity is found for an energy of ~290 keV (Sect. 13.7.21 

4.2. Long-term variability in the total emission of 
1RXSJ1 708-40 

1) In the INTEGRAL data no significant long-term time vari- 
ability for energies above 20 keV is found on time scales of 
one year and half a year in total flux and spectral index. There 
is only one measurement with a ~3cr indication for a flux in- 
crease in nine observations. The total flux (20-150keV) and 
the spectral index are stable within 22% and 18% (both lcr), 
respectivily (Sect. 1377771 MgMres[7J[2]<&'[3|i. 

2) Our analysis of a ll available ISGRI data rejects the claim 
by iGotz et al.l (120071) of correlated variability between the hard 
X-ray emission measured with ISGRI and soft X-ray emis- 
sion from 1RXS J1708-40 in relation to a (candidate) glitch 
(Sect.]3JJ\ Fig.\3$. 

4.3. Pulse profiles 

1) ISGRI measures significant pulse profiles for energies up to 
270 keV (Sect. 137277] Figures\5\&®. 

2) High-energy pulse profiles measured with XMM-Newton, 
RXTE-PCA and HEXTE and INTEGRAL-ISGRI over the 
broad energy band 0.5 ke V up to 270 keV are consistent in mor- 
phology when measured in the same differential energy bands, 
independent of the epoch of the observation; apparent differ- 
ences are either due to statistics or due to differences in energy 
response (Sect. \3.2\ Figures]6\&^. 

3) Fast morphology changes with energy are visi ble in the pulse 
profiles below ~20 keV, (earlier noted by e.g. Sugizaki et al.l 
1997tllsrael eUdl200UlKuiper et al.ll2006b In the INTEGRAL 
band the pulse profiles do not show morphology changes 
(Sect [372] Figure$6\&\7$. 

4) Three different pulse components are recognized in the pulse 
profiles: a) A hard pulse peaking around phase 0.8; which starts 
contributing to the pulse profiles above ~4 keV, b) a softer pulse 
peaking around phase 0.4; which is not apparent at hard X-rays, 
c) a very soft pulse component below 2 keV for which there is 
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an indication around phase 1 where it fills the pulse profile to a 
saw-tooth-like shape (Sect. 13.277] Fig.®. 

4.4. Total-pulsed spectra 

1) ISGRI measures the phase-averaged total-pulsed spectrum 
up to 270 keV which can be described with a power law with 
photon index T = 0.98 + 0.31 (Sect. 1337] Fig. @i. 

2) The average pulsed fraction in the INTEGRAL band is 39%. 
Note, that the pulsed flux measured above 200 keV is on the 
extrapolation of the total spectrum, consistent with being for 
100% pulsed. Below 10 keV, the pulsed fraction as measured 
with XMM-iVewfon is not constant but varies between ~3 1 % 
and ~48% in a complex way. (Sect. 13.377] Figures\4\&\9\). 

3) The total-pulsed spectrum from 0.5 keV up to 270 keV 
is very stable over time-scales of years. Firstly, the pulsed- 
emission spectra (2.7-32.1 keV) of two sets of RXTE- 
PCA data taken six years before and three years during the 
INTEGRAL operations, respectively, appear ed to be statisti - 
cally consistent within a few percent (see also lDib et alJl2008l) . 
Secondly, the summed RXTE-PCA total-pulsed spectrum joins 
smoothly that from XMM-Newton at lower energies, as well as 
those from RXTE-HEXTE and INTEGRAL at higher energies 
(Sect. UU\ F12M. 

4) Fitting INTEGRAL, RXTE-PCA and HEXTE simultane- 
ously with two power-law components yields an excellent fit 
for energies above 2.8 keV with a soft power-law model with 
photon index F = 2.79 ± 0.07 and a hard power-law model with 
photon index T = 0.86 ±0.16 (£;. =0.51, dof = 22). (Sect. 133771 
Tabled 

5) The XMM-Newton pulsed spectrum above 2.8 keV is fully 
consistent in flux and power-law index with the corresponding 
RXTE-PCA spectrum, however, there is a discontinuity in the 
XMM-Newton spectrum, a drop in flux, between 2 and 3 keV 
f Serf. 137X71 FisM. 

4.5. Phase-resolved pulsed spectra 

1) The time-averaged phase-resolved spectra, selecting three 
broad phase intervals or narrow 0.1 -phase- wide intervals, are 
connecting very smoothly over the total 0.5-300 keV energy 
band from instrument to instrument: Another confirmation of 
the long-term many-year stability of the geometry and pro- 
duction mechanisms responsible for the emission (Sect. 13.3.21 
FigM&UB- 

2) The XMM-Newton, RXTE-PCA and INTEGRAL phase- 
resolved spectra, selecting three broad phase intervals, have 
vastly different and complex spectral shapes, constituting to- 
gether the total-pulsed spectrum (Sect. 13.3.21 Fig. [§). 

3) Phlll (which contains the main hard INTEGRAL pulse) 
shows hardly any contribution to the pulsed spectrum for en- 
ergies below 2.5 keV (i.e. low 77). This part of the spectrum can 
be described by a black-body with a temperature kT of 0.23 + 
0.02 keV. Above 2.5 keV, the energy where the Phlll pulsed 
emission reaches its minimum luminosity, the spectral shape 
becomes very complex: The spectrum first hardens dramati- 
cally (and 77 increases) with a photon index F = -0.55 ± 0.14 



up to a break energy of ~ 5 keV, consistently measured with 
XMM-Newton and RXTE-PCA. After the break the spectrum 
softens to a power-law with photon index F = 1.77 ± 0.09. 
In the INTEGRAL band it hardens again to a photon index 
F = 0.86 ± 0.22 (Sect. 137372] Fig. H Tableffi. 

4) The spectra of Phi, DC-level in INTEGRAL pulse profile, 
and Phil, 'shoulder' in INTEGRAL profile and main soft-X- 
ray pulse, also appear to be complex. Both spectra are soft with 
maximum luminosity below 1 keV and exhibit discontinuities 
around 2.8 keV, the energy where we noted a drop in flux in the 
total-pulsed spectrum of XMM-Newton. The emissions seem 
to vanish above 20 keV (Sect. 13.3.21 Fig. 

5) Narrow-band phase-resolved spectroscopy in 10 phase bins 
reveals the cause of the complex behaviour of the emission 
spectra in the broader phase bins. The 10 spectra show gradual 
and sudden changes as a function of phase. The spectral com- 
plexity has become even clearer and is more extreme than seen 
in the broad-band phase-resolved spectra (Sect. 13.3.31 Fig. \ID\ . 

6) Each of the 10 pulsed spectra can be accurately described 
with a sum (free normalizations) of three components with very 
different spectral shapes: a) a soft power-law component with 
a photon index F = 3.54, b) a hard power-law component with 
a photon index F = 0.99, and c) a peculiar curved component 
which has been approximated by the sum of two logparabolic 
functions and contributes significantly between 4 and 20 keV 
(Sect. IXX31 Fig. \M Tabled. 

7) The phase distributions of the normalizations to the three 
spectral components represent three decoupled pulse profiles 
which together constitute the total pulse profile. The soft com- 
ponent shows a single peak around phase 0.4. The hard com- 
ponent shows a broad peak with a significant contribution from 
phase 0.4 to 1 . The curved component shows a narrower peak 
(~0.25 in phase FWHM) which has its maximum in phase in- 
terval 0.7-0.8 (Sect. 1X3731 Fig.TB. 

5. Discussion 

The detailed analysis performed in this work on the high- 
energy emission (0.5-300 keV) from 1RXS J1708-40 led 
to similar conclusions as were reached earlier for AXP 



4U 0142+61 dden Hartog et al.ll2008l) . Namely, that genuinely 
different pulse components with different spectra contribute to 
the total-pulsed emission. In the case of 1RXS J1708-40 not 
only the variation with energy of the structure of the pulse pro- 
files could be used to identify these different components, sim- 
ilar to the case of 4U 0142+61, but also the apparent fine struc- 
ture in the phase-resolved spectra provided an efficient means 
to dise ntangle the pulse d emiss ion into different components. 

In |d en Hartog et al .1 (120081) we discussed the implications 
of the derived high-energy characteristics of 4U 0142+61 for 
models developed in recent years for explaining the non- 
thermal, luminous hard X-ray emission from AXPs: 1) a 
quantum electro-dynamic model (OE D) by iHev l & Hernquist 
(200533), 2) a corona model by Beloborodov & Thompson 
d2007l) . further elaborated bv lLvubarskv & Eichlerl (120081) an d 
3) a resonant upscattering model by Barin g~ Harding! (12007 ). 
We concluded that all three models have difficulties in explain- 
ing many features of the detailed results (See the discussion in 
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the paper bv lden Hartog et al]|2008l) . A confrontation of these 
models with our new results for 1RXS J1708-40 leads to sim- 
ilar conclusions. It is apparent that any attempt to explain the 
high-energy emission (0.5-300 keV) should consider a stable 
(at least over a decade) three-dimensional geometry taking all 
relevant angles (the angle between the magnetic and the spin 
axis and the viewing angle) into account, as well as the physical 
production processes taking place on the surface of the neutron 
star and at different sites in the atmosphere and magnetosphere. 

For 1RXS J1708-40 our results reveal a clear separation in 
phase, and therefore also of the production sites, of all pulsed 
hard X-ray emission above 20 ke V with photon index F ~ 1 .0 
from the pulsed soft thermal emission below 4 keV. This soft 
emission can be described with a black-body spectrum with 
kT < 0.4 keV and a power-law spectral tail with photon in- 
dex T ~ 3.5. Such a clear separa tion i s consi stent with pre- 
dictions from the QED model by iHeyll ( 2007 ). However, the 
observed total and the pulsed hard X-ray spectra are harder 
(index ~1.0) than the photon index of 1.5 predicted by this 
model. Contrarily, the clear separation in phase between the 
thermal and non-thermal pulses makes an interpretation of the 
hard X-ray emission in terms of resonant magnetic upscat- 
tering of the thermal soft photons by ultra-relativistic elec- 
trons, accelerated along either open or closed magnetic field 
lin es unlikely (proposed in the Compton-upscattering model 
by Baring & Harding! 2007 ). At least part of the thermal tar- 
get photons should also arrive in the narrow phase interval of 
the non-thermal pulse. This problem might be circumvented if 
the thermal target photons are part of the isotropic thermal DC 
emission, and the thermal pulse at phase 0.4 is not related to 
the hard X-ray emission. 



As mentioned before, iGotz et alj (12007b reported the de- 
tection of correlated variability in hard X-ray flux and spectral 
index with similar correlated variations at low-energy X-rays 
from 1RXS J1708-40. This would imply correlated production 
processes at the same sites in the magnetosphere. However, 
our analysis disproved the claimed variability in the ISGRI- 
detected hard X-ray flux. 

Very peculiar is the spectral shape of the component which 
contributes significantly to the total-pulsed emission from 
1RXS J 1708-40 between 4 keV and 20 keV. It is not obvi- 
ous what underlying emission processes can produce such a 
spectral shape. This component exhibits the most narrow pulse 
profile with a width of ~0.25 in pulsar phase (FWHM, see 
Fig. [TTl l. It appears to be completely separated from the soft 
pulse, but it is centered in phase on the hard X-ray pulse which 
is about twice as broad. This means that the production site 
in the magnetosphere of this curved component is confined 
to a narrow region, possibly directed along open or the last 
closed magnetic-field lines, sim ilar to different dipole scenarios 
proposed for radio pulsa r s (e.g.lRuderman & Sutherland] 19751: 



Usov & Melrose 1995; Harding & Muslimo 
1983 1 ; iMuslimov & Harding|l2003t ICheng et 



ov 
al. 



1998 



Arons 



Romani 



1996; Hirotani 2006). Like is the case for radio pulsars, this ge- 
ometry appears to be very stable. The RXTE monitoring obser- 
vations proved that the pulse profiles and phase resolved spec- 
tra were stable over 9 years. This might indicate that the non- 
thermal emission is produced higher-up in the magnetosphere, 



above the region closer to the neutron star with a strong twisted 
(toroidal) field. 

The detailed results we obtained for AXPs 4U 0142+61 
and 1RXS J1708-40 clearly prove that the total persistent high- 
energy emission from AXPs consists of DC emission and 
pulsed emission ex hibiting different spectra (already shown in 
iKuiper et al.l l2006). Furthermore, the pulsed emission consists 
of multiple components with different spectra. As a conse- 
quence it does not seem meaningful to make detailed model fits 
to just the total AXP spectra to derive a physical interpretation 
of the emission. The physical origin and production site of each 
different component has to be considered. Examples of model 
fitting to the total X-ray spectra of magnetars are recent ap- 
plication s of the resonant cyclot ron scattering model first pro- 
posed by Thompson et alj J2002f). Recent developments of this 



model can be found in Lyutikov & Gavriil (2006), Giiver et al 



(120061 120071) and lFernandez & Thompson! (120071) . This model 
computes the result of multiple scatterings of soft photons in 
the magnetospheres and was developed to explain the appar- 
ently thermal black-body spectra with a power-law extension 
to higher energies, as me asured for energie s below ~10 keV. 
In the extensive work by iRea et al . (2008) detailed fits with 
this model are merely made to the high-energy total spectra 
of a comprehensive set of magnetars, including 1RXS J1708- 
40, leading to interesting estimates of e.g. optical depths and 
magnetospheric electron densities at the resonances. However, 
our temporal and spectral decoupling of different components 
in the pulsed emission (Figures [TOl & fTTb. shows that between 
~4 keV and 20 keV a peculiar curved spectral component con- 
tributes significantly to the pulsed emission, and is unrelated 
to the soft thermal emission. This component contributes also 
significantly to the total emission detected from 1RXS J1708- 
40 in the energy interval where the resonant upscattering is 
effective. Hence, it will obviously have an adverse impact on 
the parameter estimations if this component has not been sub- 
tracted. On the other hand, the resonant cyclotron-scattering 
model might offer a good explanation for the soft pulse peaking 
at phase 0.4. Spectral modeling of the corresponding spectrum, 
and separately of the thermal spectrum of the DC emission will 
probe the model parameters in different sites in the magneto- 
sphere. 

An important parameter to be reproduced by the models is 
the photon index of the power-law spectral shape in the hard 
X-ray window and the cut-off energy. The index F ~ 1 .0 of the 
non-thermal pulse and a spectral break at a few hundred keV is 
consistent with a bremsstrahlung origin as first prop osed in the 
corona model by Beloborodov & Thompson! ( 2007b. But, fur- 
ther e laboration of the corona model by iLvubarskv & Eichlei 
( 2008) lead to a different hard X-ray spectrum, namely a super- 
position of bremsstrahlung radiation from the hot atmosphere 
and Comptonization radiation from the extended corona. The 
latter produces a hard spectrum up to the MeV band and is 
inconsistent with the COMPTEL upper limits at MeV ener- 
gies. Furthermore, the stability of the configuration over 9 years 
is a challenge f or the corona mode l . For the resonant upscat- 
tering model of Baring & Harding! (12007b the predicted pho- 
ton indexes are considerably softer and the spectra extend into 
the MeV range, possibly to significantly higher gamma-ray 
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energies. In this respect we should emphasize that we donot 
see any indication for a break in the INTEGRAL spectrum 
(<300 keV) of 1RXS J1708-40, more importantly, we donot 
know for 1RXS J1708-40 and 4U 0142+61 what the shapes of 
the spectra are above a few hundred keV. Extreme cases are al- 
lowed by the data, e.g. a hard spectrum with index F ~ 1 .0 with 
a sharp cut-off, as well as with a break or bend with a flat spec- 
trum extending underneath the COMPTEL upper limits into the 
MeV or even GeV regions. 

Similarly hard spectra with F ~ 1.1 have been measured 
in the hard X-ray window up to 60-200 keV for several young 
rotation-powered pulsars, e.g. AX J1838.0-0655, PSR J1846- 
0258, PSR J181 1-1925, PSR J1930+1852 and PSR J1617- 



5055 dKuiperetal.1 120081: iKuiper & Hermsenll2008l) . We also 
know that these hard power-law spectra donot extrapolate to the 
MeV range (no detections by COMPTEL). Furthermore, those 
young pulsars detected with COMPTEL in the MeV range, 
Crab (for a multi-w avelength analysis simi lar to this work se e 
Kuiperet allboOlh and PSR B1509-58 dKuiperet all 1 19991) . 
have flat (F ~ 2) spectra at MeV energies, and started with 
harder spectra at hard X-ray energies. In addition, for three mil- 
lisecond pulsars also an index F ~ 1.1 has been measured in 
the X-ray band up to 25 keV (PSR J0218+4232, PSR B1821- 



24 an d PSR B1937+21: IKuiper et al.ll2004bl: ICusumano et al 



20031) . For one of these millisecond pulsars, PSR J0218+4232, 
detection has been claimed even be tween 300 MeV and 1 GeV 
with a soft spectral index of ~2.6 dKuiper et al. 2000), thus a 
drastic spectral bend or break is required between the hard X- 
rays and the high-energy gamma rays. These results on young 
as well as on old millisecond radio pulsars prove that in dipole 
fields such hard non-thermal spectra can be produced in stable 
configurations. The major difference with the AXPs is that the 
energy source for ordinary pulsars is rotational energy and for 
AXPs/magnetars magnetic energy. It is interesting to note that 
recently during a ~50-day outburst magnetar-like bursts have 
been reported from a very young rotation-p owered pulsar ( PSR 
J 1846-025 8) with a stron g magnetic field dKumar & Safi-Harbl 
2008HGavriil et al.ll2008l) . 

A possible scenario for the persistent high-energy emis- 
sion from AXPs seems to emerge in which a soft thermal 
component, for 1RXS J1708-40 peaking at phase 0.4, origi- 
nates in regions above and close to the polar cap of the neu- 
tron star with resonant cyclotron scattering offering a physi- 
cal explanation. A large fraction of the thermal photons will 
have a more isotropic distribution. These photons can be re- 
sponsible for the soft part of the DC emission spectrum, also 
with a resonant cyclotron-scattering origin, but with e.g. dif- 
ferent optical depths and magnetospheric electron densities at 
the resonances. For future studies, it will be interesting to ac- 
curately derive a spectrum of just the DC emission for model 
fitting. Furthermore, the non-thermal emission might be pro- 
duced higher up in the outer magnetosphere, where the mag- 
netic field strengths is in the range estimated for normal radio 
pulsars. The required strong current of electrons and positrons 
can be sustained by the outward transfer of magnetic helicity 
from the inner zone where it is injected. Then the accelera- 
tion of particles and the production of the observed non-thermal 
photons has to occur in stable narrow zones along the curved 



field lines to produce e.g. the narrow hard X-ray pulse seen 
for 1RXS J1708-40. In the outer magnetosphere the scenario 
might resemble those proposed for the high-energy emission 
from rotation-pow e red p ulsars. In fact, the approach taken by 
iBaring & Hardingl d2007l) is a first attempt along these lines. 
However, they concluded that the resonant magnetic Compton 
upscattering is most effective close to the magnetar surface. 

From only two of the 13 confirmed magnetars radio emis- 
sion has been detected and their radio emission properties 
have been measu red in great det a il rece n tly, notably from AX P 
XTE Jl 8 10- 197 lHalpern et al.1 d2005h : ICamilo et al.1 d2006l): 
Kramer et al.1 (|2007|) and AXP IE 1547.0-5408 dCamilo et al 



20071 12008). These very interesting AXPs show transient ra- 
dio emission , which is apparentl y related to the X-ray vari- 
it f 



ability (e.g. Ilbrahim et all 12004 iGotthelf & Halpernl 12005 
iGelfand & Gaenslerl 120071 lHalpern et al.1 12008b . These radio 



AXPs appear to be very highly linearly polarized and have very 
flat spectra over a wide range of radio frequencies, very similar 
characteristics to th ose of many young ordinary radio pulsars. 



Camilo et al.l d2008l) note that, while unexpected a priori, it ap- 



pears that ideas developed to understand the geometry of emis- 
sion from ordinary pulsars along open dipolar field lines also 
apply to known radio magnetars. However, they also conclude 
that both radio AXPs have also variable pulse profiles and radio 
flux densities, not typical for young pulsars. 

Upon submission of this work, an extensive treatment 
appeared of pair creation processes in an ultra-strong mag- 
netic field and par ticle heating in a dynamic magnetosphere 
dThompsonl 120 08 a), as well as a self-consistent model of the 
inner accelerator w ith implications for pulsed radio emission 
(Thompson 2008b). In the first paper Thompson considers the 
details of QED processes that create electron-positron pairs in 
magnetar fields. He discusses how the pair creation rate in the 
open-field circuit and the outer magnetosphere can be strongly 
enhanced by instabilities near the light cylinder and how this 
very high rate of pair creation on the open magnetic field lines 
can help to stabilize a static twist in the closed magnetosphere 
and to regulate the loss of magnetic helicity by reconnection at 
the light cylinder. In the second paper, one of the possible sce- 
narios for explaining the rapid radio variability, broad pulses 
and hard radio spectra comprises a strongly variable current 
in the outer magnetosphere and a high rate of pair creation sus- 
tained by a turbulent cascade. It is interesting to note that stable 
geometries are being considered to explain the production of 
transient radio emission for AXPs involving high rates of pair 
creation on the open magnetic-field lines in the outer magne- 
tosphere. Such scenarios might also be required to explain the 
production in the outer magnetosphere on the orders of magni- 
tude more luminous non-thermal X-ray emission. 
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